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In the framework of phenomenological static and QCD-motivated model potentials 



' for heavy quarkonium, we compute the bc mass spectrum as well as its IS hyperfine 

splitting using the recently introduced shifted I- expansion technique. We also predict 
^ ■ the leptonic constant fs c of the lightest pseudoscalar B c and fs* of the vector B* 

in 

■ states taking into account the one-loop and two-loop QCD corrections. Further, 

we use the scaling relation to predict the leptonic constant of the nS-states of the 



bc system. For the sake of comparison, we use the same fitting parameters of our 
previous model potentials. From our results we conclude that shifted /-expansion 



^ ■ method has the same accuracy, convergence and status as our previous work. 
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1. INTRODUCTION 

The spectrum and properties of the be quarkonium system have been calculated various 
times in the past in the framework of heavy quarkonium theory [1]. Moreover, the recent 
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discovery [2] of B c meson (the lowest pseudoscalar 1 S state of the B c system) opens up 
new theoretical interest in the subject [1,3-8]. The Collider Detector at Fermilab (CDF) 
Collaboration quotes M Bc = 6A0%%f 3 GeV [2]. This state should be one of a number of 
states lying below the threshold for emission of B and D mesons. Further, such states are 
very stable in comparison with their counterparts in charmonium (cc) and bottomonium (bb) 
systems. A particularly interesting quantity should be the hyperfine splitting that as for cc 
case seems to be sensitive to relativistic and subleading corrections in the strong coupling 
constant a s . For the above reasons it seems worthwhile to give a detailed account of the 
Schrodinger energies for cc, bb and be meson systems below the continuum threshold. 

Because of the success of the nonrelativistic potential model and the flavour independence 
of the q ± q 2 potential, we choose a wide class of phenomenological and a QCD-motivated po- 
tentials by insisting upon strict flavor-independence of its parameters. We also use a poten- 
tial model that includes running coupling constant effects in both the spherically symmetric 
potential and the spin- dependent potentials to give a simultaneous account of the properties 
of the cc, bb and be heavy quarkonium systems. Since one would expect the average val- 
ues of the momentum transfer in the various quark-antiquark states to be different, some 
variation in the values of the strong coupling constant and the normalization scale in the 
spin-dependent potential should be expected. To minimize the role of flavor-dependence, we 
use consistent values for the coupling constant and universal QCD renormalization scale for 
each of the levels in a given system. 

Recently, Kwong and Rosner [7] predicted the masses of the lowest vector and pseu- 
doscalar states of the be system using an empirical mass formula and a logarithmic poten- 
tial. Eichten and Quigg [1] gave a more comprehensive account of the energies and decays 
of the B c system that was based on the QCD-motivated potential of Buchmuller and Tye 
[9]. Gershtein et al. [8] also computed the energies and decays of the B c system using the 
QCD sum-rule calculations. Baldicchi and Prosperi [6] have fitted the entire light-heavy 
quarkonium spectrum and computed the be spectrum based on an effective mass operator 
with full relativistic kinematics. Fulcher [4] extended the treatment of the spin-dependent 
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potentials to the full radiative one-loop level using the renormalization scheme developed 
by Gupta and Radford [10]. Ebert et al. [1] investigated the B c meson masses and decays 
in the relativistic quark model. Very recently, we reproduced the cc, bb and be spectroscopy 
by applying the shifted large-N expansion technique (SLNET) on the nonrelativistic and 
relativistic wave equations using a group of static and improved QCD motivated potentials 
[11]. 

Encouraged by the success of SLNET application on heavy quarkonium systems [11,12], 
we extend the previous works in [11] by applying the shifted /- expansion technique (SLET) 
[13] on the Schrodinger equation to reproduce the be spectroscopy. The SLET has been re- 
cently introduced to solve mathematically the two and three-dimensional Schrodinger equa- 
tion for spherically and cylindrically symmetric potentials [13]. Further, the KG and Dirac 
equations with radially symmetric Lorentz vector and Lorentz scalar potentials have also 
been solved [13]. We anticipate that, by working SLET, we would be able to obtain a better 
understanding on the status, convergence and accuracy of this method among the other 
methods [11]. We also consider it as a complementary investigation to our previous works 
[11]. Further, our results will enable us to check clearly that SLET is just a parallel per- 
turbative expansion method, i.e., /-expansion procedure in a similar manner to /c-expansion 
procedure. 

The contents of this article are as follow: in section 2, we present the solution of the 
Schrodinger equation using the SLET for the the non-self conjugate Qq meson mass spec- 
trum. In section 3 we briefly present the potentials used. In section 4 we give the first-loop 
and second- loop correction of the B c leptonic decay constant. Finally, a discussion and 
conclusion appear in section 5. 

2. SCHRODINGER WAVE EQUATION 

In previous papers [11,12] we have applied the shifted 1/N expansion technique (SLNET) 
to solve nonrelativistic and relativistic wave equations. The method starts by writting 
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the original wave equation in an N- dimensional space which is sufficiently large and using 
expansion 1/k as a perturbation parameter [14]. Here k = N + 21 — a, N being the number 
of spatial dimensions of interest, / the angular quantum number, and a is a suitable shift as 
an additional degree of freedom and is responsible for speeding up the convergence of the 
resulting energy series. In this work another method called the shifted /-expansion technique 
(SLET) which is simply consists of using 1/1 as an expansion parameter, where I — I — a, 
I is an angular quantum number, a is a suitable shift which is mainly introduced to avoid 
the trivial case / = 0. The choice of a is physically motivated so that the next to the leading 
energy eigenvalue series vanish as in SLNET. The method does not require writting the 
N- dimensional form of the wave equation and we expand it directly through the quantum 
number / involved in the problem. This method seems more flexible and simpler in treatment 
and has a quite different mathematical expansion than SLNET. Like SLNET, SLET is also a 
pseudoperturbative technique. The radial part of the Schrodinger equation for an arbitrary 
spherically symmetric potential V(r) (in units H — 1) 

f 1 d 2 1(1 + 1) Tr , A , . „ , . 

where \i = (m 9 mg) /(m q + mq) is the reduced mass for the two constituent interacting 
particles and E n £ denotes the Schrodinger binding energy. Equation (1) can be rewritten as 



1 d 2 

+ 



f + (2a + 1) I + a (a + V 



+ V(r)\u(r) = E ntl u(r), (2) 



4/x dr 2 4/ir 2 

where I — I — a with a representing a proper shift to be calculated later. It is clear from Eq. 
(2) and our previous work [11] [cf. Eq. (1) in hep-ph/0303182] that I is as much as k/2 of 
Imbo et al. [14]. We follow the shifted /-expansion method [13] by defining 

and also the energy eigenvalue expansion [11,12] 

oo 7( 2 - m ) 

E n ,i =E -n~ Em - (4) 

m=0 



-1/2 

Here y = I (r/r — 1) with r is an arbitrary point where the Taylor expansions is being 
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performed about and Q is a scale to be set equal to I at the end of these calculations. 
Inserting Eqs. (3) and (4) into Eq. (2) yields 

(-l) m (m + l)y m 



Id 2 1 / 7 ,„ a(a + l) 

+ — [l + (2a + 1) + 



4/i dy 2 4/x 



I 



E 

m=0 



jm/2 



,2 oo 



'd m \/(r )\ (r y) m -(2-m)/2 



m=0 



m! 



(y) = UrXn r (y)- 



(5) 



Hence the final analytic expression for the l/l expansion of the energy eigenvalues appro- 
priate to the Schrodinger particle is [13] 



r 2 00 ci ^ 
' n x y(l— m) 



^ m=0 



-Em. 



(6) 



Now we formulate the SLET for the nonrelativistic motion of spinless particle bound in 
spherically symmetric potential V(r). Equation (6) is exactly of the type of Schrodinger 
like equation for the one dimensional anharmonic-oscillator and has been investigated for 
spherically symmetric potential by Imbo et al. [14]: 



1 r 2 V(r ) 
4fi Q 



+ 



a(a + 1) 
4/i 



+ 7 



1. (2a + 1) 



(7) 



where the expressions 7^ and 7^are given explicitly in Appendix A. Thus, comparing Eq. 
(6) with Eq. (7) gives 

Q 



E = V(r ) + 



4 / ur 2 , 



2 ' 



(8) 



Ei — -p 



1\ (2a + 1) 
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and 
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The quantity r is chosen so as to minimize the leading term, E , that is, [12] 

— — = and > 0, 

ar (Ltq 

which yields the relation 
Q = 2/ir V / (r ). 



(12) 



(13) 



Further, to solve for the shifting parameter a, the next contribution to the energy eigenvalues 
is chosen to vanish [11-14], i.e., E x = 0, which provides smaller contributions for the higher- 
order corrections in (4) compared to the leading term contribution (8). It implies that the 
energy states are being calculated by considering only the leading term E , the second-order, 
E 2 and the third-order, E 3 corrections. Therefore, the shifting parameter is determined via 
[1 + 2/i(2n r + 



a = 



(14) 



with 



1 

2/1 



3 + 



r V"(r ) 
V'(r ) 



1/2 



(15) 



Thus, the Schrodinger binding energy (4) to the third order is 



E n ,i = V(r ) + - 

Further, setting I - 
satisfies r as 



4/i I \l\ 



(16) 



Q which rescales the potential, we derive an analytic expression that 



21 + 



1 + (2n r + 1) 



3 + 



r V"(r ) 
V'(r ) 



1/2 



[tyr*V'{r ) 



1/2 



(17) 



where n r = n — 1 is the radial quantum number. Once r is being found through Eq. 
(17) for any arbitrary state, the determination of the binding energy for the Qq system 
becomes relatively simple and straightforward. Finally, the Schrodinger binding mass can 
be determined via 
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M(Qq) =m q + m Q + 2E n +. (18) 

It is being found that for a fixed n, the computed energies become more accurate as I 
increases [11-14]. This is expected since the expansion parameter 1/1 becomes smaller as I 
becomes larger since the parameter 7 is proportional to n and appears in the denominator 
in higher-order correction. 

3. SOME MODEL POTENTIALS 

The be system that we investigate is often considered as nonrelativistic system [1] and 
consequently our treatment is based upon Schrodinger equation with a Hamiltonian [1,4] 

H = -^- + V(r) + V SD: (19) 

where we have supplemented our nonrelativistic Hamiltonian with the standard spin- 
dependent terms [1,11,15] 

*° - v - = S<-^ 3 M- < 20 > 

Here, the spin dependent potential is simply a spin-spin part [1,15] that would enable us 
to make some preliminary calculations of the energies of the lowest two S-states of the be 
system. The potential parameters in this section are all strictly flavor-independent and 
fitted to the low-lying energy levels of cc and bb systems. Like most authors (cf. [1]), we 
determine the coupling constant a^m;?) 1 from the well measured hyperfine splitting for the 
lS(cc) state [16] 

A£ HF (1S, exp) = M m - M Vc = 117.2 ± 1.5 MeV, (21) 
and for the 2S(cc) state [16-18] 



Due to the lack of any experimental splitting data on the B c meson, as there is one established 
state, we have fitted the coupling constant to reproduce the available cc splittings. 
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A£ HF (2S, exp) = My - M Vc , = 32 ± 14 MeV, (22) 

for each desired potential to produce the center-of-gravity (cog) of the M^(1S) value. The 
numerical value of a s is found to be dependent on the potential form and also be compatible 
with the other measurements [1,3,4,6-8]. Therefore, the lS-state hyperfme splitting [1,11,15] 
is given by 2 

AE UF = -^ \R 19 (0) | 2 , (23) 
9m c m b 

with the radial wave function at the origin is determined via [15] 

|iM0)|* = 2„(^). (24) 
Hence, the total mass of the low-lying pseudoscalar B c meson is [11] 

M Bc (0-) =m c + m b + 2E lfi - 3 AE UF /4, (25) 
and for the vector B* meson 

M B ,(1-) = m c + m b + 2E lfi + AE UF /A. (26) 
Hence, the square-mass difference can be simply found as 

AM 2 = Mf.(l-) - M 2 Bc {0~) = 2AE UF [m c + m b + 2E lfi - AE RF /4] . (27) 

The perturbative part of such a quantity was evaluated at the lowest order in a s . Baldicchi 
and Prosperi [6] used the standard running QCD coupling expression 

Atx 

as{Q) = Ai 2 N, ( Q 2 \ ' ^ 
(ll-fn/J ln(^) 



2 To the moment, the only measured splitting of nS— levels is that of rj c and J/ip, which allows us 
to evaluate the so-called SAD using M^(1S) = (3Mj/^, + M?7 c )/4 and also M(nS) = My (nS) - 



{M m - M Vc )/4n [15,19]. 
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with n/ = 4 and A = 0.2 GeV cut at a maximum value a s (0) = 0.35, to give the right 
J/ij} — r) c splitting (21) and to treat properly the infrared region [6]. Further, Brambilla 
and Vairo [3] took in their perturbative analysis 0.26 < a s (/i = 2 GeV) < 0.30. After 
the observation of the r/ c (2S) meson [17], Badalian and Bakker [18] determination of the 
coupling constant aup(^i, IS) ^ 0.335 is rather large with ^ ~ ^M{J/ip) = 1.55 GeV and 
ttHF(/^2, 2S) ~ 0.18 with ji 2 — 2M(if)') = 7.4 GeV which implies the very large value for the 
renormalization scale. They [20] also used a s (/j, = 0.92 GeV) ~ 0.36 for all states, but the 
splittings do practically not change if a s (/j, = 1.48 GeV) = 0.30 is taken. This result shows 
that the strong coupling constant depends on the renormalization scale and it changes from 
one state to another [18]. Further, Motyka and Zalewski [21] found a s (m 2 c ) = 0.3376 and 
from which they calculated a s (ml) = 0.2064 and a s (4/i^ c ) = 0.2742. Therefore, it is quite 
clear that the coupling constant is dependent on the quarkonium system. 

The commonly used potentials are of two types: (i) pure phenomenological and (ii) partly 
phenomenological, but motivated by both perturbative and non-perturbative QCD at short 
and long distances. In this work we use the following types of potentials: 

A. Static potentials 

It is seen that the most phenomenological potentials in Eq. (19) may be gathered up in 
general form [21,22]: 

V(r) = -ar~ a + br p + c 0<a,/3<l, a,b>0. (29) 

where c may be of either sign. The mixed powerlaw (29) comprises more than ten potentials 
given by Refs. [21,22]. It seems quite reasonable that the short range behaviour of the 
quarkonium potential is less singular than — r _1 and the confining potential does not rise 
quickly as r due to the screening effects of quark pair creation. Therefore, the effective 
quarkonium potential consist of two terms, one of which, V v (r) = —ar~ a , transforms like a 
time-component of a Lorentz 4- vector and the other, V s (r) = br 13 + c, like a Lorentz scalar. 
We limit our study to the a — ft — v case 
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V(r) = 



-ar~ v + br u + c < v < 1, a,6>0. 



(30) 



proposed by Lichtenberg et al. [23]. Like most authors [1,4,11,23], we consider a class of 
static potential which give reasonable accounts for the cc and bb spectra. This comprises 
a wide class of potentials presented explicitly in our previous works [11] some are QCD- 
motivated potentials like Cornell {v — 1), and other typical pure phenomenological potentials 
like Song-Lin {y = 1/2), Turin [y = 3/4), Logarithmic [y — > 0) which are all belonging to 
the class (30) and Martin (a — 0, (5 — 0.1) belonging to the class (29). The motivation of 
this choice is that all of these potentials lie very close together in the range of distances 
0.1 < r < 1, which is the characteristic interval of cc and bb spectra. 



B. QCD-motivated potentials 

1. Igi-Ono potential 

Buchmuller and Tye [9] proposed a potential which is consisting of two parts, at short 
distances the two- loop perturbative calculation of the interquark one-gloun exchange [24]: 

16tt 1 



v OGE V 



with 



f(r) = In 



25 rf{r) 



1 462 lnf(r) | 2 lE + § ^ 



625 f(r) 



f(r) 



(31) 



r 2 A- — 

MS 



+ b 



(32) 



where rif = 4 is the number of flavors with mass below /i and 7^ = 0.5772 is the Euler's 
number. Moreover, at long distances the interquark potential grows linearly leading to 
confinement as 



V L {r) = ar. 



(33) 



Therefore, the Igi-Ono potential is [24] 



\/W= 4 )( r ) = V^ i] + ar + dre- 9r , 



(34) 
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where the term dre~ gr in (34) is added to interpolate smoothly between the two parts and 
to adjust the intermediate range behavior by which the range of Ajjg is extended keeping 
linearly rising confining potential. Numerical calculations show that potential is good for 
A^r in the range 100-500 MeV keeping a good fit to the cc and bb spectra. The QCD 
coupling constant a s in (20) is defined in the Gupta-Radford (GR) renormalization scheme 
[10] 

a s = t r-, (35) 

(33-2„ / )ln( x H_) 

where the scale parameter in the Gupta-Radford (GR) renormalization scheme Aqr [10] is 



related to Aj^ by 



A GR = A MS ex P 



49 - 10ra//3 



_2(33-2n/) 

Thereby, the three types of this potential are displayed in [11]. 



(36) 



2. Improved Chen-Kuang potential 



Chen and Kuang [25] proposed two improved potential models so that the parameters therein 
all vary explicitly with Ajjg so that these parameters can only be given numerically for 
several values of Aj^. Such potentials have the natural QCD interpretation and explicit A^r 
dependence both for giving clear link between QCD and experiment and for convenience in 
practical calculation for a given value of Aj^. It has the general form 



\/( n / =4 )( 



kr 



16vr 1 
25 rf(r) 



462 In f(r) , 2 lE + f 



1 - 



625 f(r) 



f(r) 



(37) 



where the string tension is related to Regge slope by k = The function f(r) in (37) is 

1 2 



f(r) = In 
with 



A 



+ 5.10- A(r) 



MS' 



(38) 



A(r) 



1 - 



1A 



MS 



4Ai— 

MS. 













1 — exp | — 
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A 



(39) 



MS' 
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The scale parameter A-^ is very close to the value of determined from the two-photon 
processes and is also close to the world-averaged value of A-p^. The fitted values of its 
parameters are displayed in Ref. [11]. 



4. LEPTONIC CONSTANT OF THE 5 c -MESON 



The study of the heavy quarkonium system has played a vital role in the development of 
the QCD. Some of the earliest applications of perturbative QCD were calculations of the 
decay rates of charmonium [26]. These calculations were based on the assumption that, in 
the nonrelativistic (NR) limit, the decay rate factors into a short-distance (SD) perturbative 
part associated with the annihilation of the heavy quark and antiquark and a long-distance 
(LD) part associated with the quarkonium wavefunction. Calculations of the annihilation 
decay rates of heavy quarkonium have recently been placed on a solid theoretical foundation 
by Bodwin et al. [27]. Using NRQCD [28] to seperate the SD and LD effects, Bodwin et al. 
derived a general factorization formula for the inclusive annihilation decay rates of heavy 
quarkonium. The SD factors in the factorization formula can be calculated using pQCD [19], 
and the LD factors are defined rigorously in terms of the matrix elements of NRQCD that 
can be estimated using lattice calculations [5]. It applies equally well to S-wave, P-wave, 
and higher orbital-angular-momentum states, and it can be used to incorporate relativistic 
corrections to the decay rates. 

In the NRQCD [28] approximation for the heavy quarks, the calculation of the leptonic 
decay constant for the heavy quarkonium with the two-loop accuracy requires the matching 
of NRQCD currents with corresponding full-QCD axial-vector currents 



^ NRQCD = -Xl*** J " QCD = *»V, («) 

where b and c are the relativistic bottom and charm fields, respectively, xl an d V'c are the 
NR spinors of anti-bottom and charm and v A is the four-velocity of heavy quarkonium. The 
NRQCD lagrangian describing the £? c -meson bound state dynamics is [30] 
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£nrqcd = Aight + 4 (iD + D 2 /(2m c )) ^ + X \ (*A> - D 2 /(2m 6 )) Xb + • • • , (41) 

where £u g ht is the relativistic lagrangian for gluons and light quarks. The two-component 
spinor field ip c annihilates charm quarks, while Xb creates bottom anti-quarks. The relative 
velocity v of heavy quarks inside the 5 c -meson provides a small parameter that can be used 
as a nonperturbative expansion parameter. To express the decay constant f Bc in terms of 
NRQCD matix elements we express J A in terms of NRQCD fields ip c and Xb- The A = 



QCD 



current-component contributes to the matrix element and consequently the J x 
following operator expansion 



has the 

QCD 



<0|&7 A 75C|£ C (P)> =if Bc P\ (42) 

where \B C (P)) is the state of the 5 c -meson with four-momentum P. Only the A = compo- 
nent contributes to the matrix element (42) in the rest frame of the S c -meson. It has the 
standard covariant normalization 

7^3 / r Bc (p')^B c (p)d 3 p = 2E5^(p> - p), (43) 

and its phase has been chosen so that f Bc is real and positive. Hence, the matching yields 

6 7 75 c = K xliP c + K 2 (B Xb y.-Dij c + • • • , (44) 

where K = K Q (m c ,mb) and K 2 = K 2 (m c ,mi ) ) are Wilson SD coefficients. They can be 
determined by matching perturbative calculations of the matrix element (0| &7°7 5 c \B C ) , a 
contribution is mostly coming up from the first term in 



<0| 6 7 °7 5 c \B C ) ^ = K (0| xUc \B C ) + K 2 (0| (D Xb )t.D^ c \B, 



QCD 



+•••, (45) 

NRQCD 



NRQCD 

where the matrix element on the left side of (45) is taken between the vacuum and the state 
\B C ) . Hence, equation (45) can be estimated as: 

(0\ X t^ c \B c ) 2 ^^\R 1S (0)\ 2 . (46) 

7T 

Onishchenko and Veretin [30] calculated the matrix elements on both sides of equation (45) 
up to a 2 s order. In one-loop calculation, the SD-coefficients are 
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Kr> — 1 and K? = 



8/i 2 ' 



(47) 



with /i denned after Eq. (1). Further, Braaten and Fleming in their work [31] calculated 
the perturbation correction to K up to order- a s (one-loop correction) as 



1 + Ci- 



7T 



(48) 



with ci being calculated in Ref. [31] as 
2 - 



ci = - 



m b -m c . m b 
in — 



m 6 + m c 

Finally, the leptonic decay constant for the one-loop calculations is 



(49) 



(l-loop) 



1 - 



a s (fi) 



7Y 



m b - m c m b 
2 in — 



m b + m c 



NR 



(50) 



where the NR leptonic constant [32] is 
#is(0)| 



NR 



TTMf 



(51) 



and /x is any scale of order m b or m c of the running coupling constant. On the other hand, 
the calculations of two-loop correction in the case of vector current and equal quark masses 
was done in Ref. [33]. Further, Onishchenko and Veretin [30] extended the work of Ref. 
[33] into the the non-equal mass case. They found an expression for the two-loop QCD 
corrections to £? c -meson leptonic constant given by 



K (a s , M/fi) = 1 + Cl (M//i)^^J + c 2 (M//i) 

7T \ 71 



+ 



(52) 



where ci{Mj ii) is the two- loop matching coefficient and with ci j2 are explicitly given in Eq. 
(49) and (cf. Ref. [30]; Eqs. (16)-(20) therein), respectively. In the case of £? c -meson and 
pole quark masses (m b = 4.8 GeV, m c = 1.65 GeV), they found 



(2-loop) 
B c 



1 - 1.48 



'a s (m b y 



- 24.24 



'a s (m b ) 



7Y 



7Y 



NR 
B c ■ 



(53) 



Therefore, the two-loop corrections are large and constitute nearly 100% of one-loop correc- 
tion as stated in Ref. [30]. where /i is chosen as the scale of the running coupling constant. 
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5. DISCUSSION AND CONCLUSION 



We use Eq. (21) to determine the position of the charmonium center-of-gravity M^,(1S) mass 
spectrum. Further, we fix the coupling constant a s (m c ) for each potential. For simplicity we 
neglect the variation of a s with momentum in (28) to have a common spectra for all states 
and scale the splitting of be and bb from the charmonium value in (21). The consideration of 
the variation of the effective Coulomb interaction constant becomes especially essential for 
the T particle, for which a s (T) ^ a s (ip). 3 We follow our previous works [11] to calculate the 
corresponding low-lying center-of-gravity M T (1S) and consequently the low-lying M Bc (lS). 

Table 1 reports our prediction for the Schrodinger mass spectrum of the four lowest cb 
S-states together with the first three P- and D-states below their strong decay threshold 
for different static potentials. Since the model is spin independent and as the energies of 
some singlet states of quarkonium families have not been measured [11,19,23], a theoretical 
estimates of these unknown levels introduces uncertainty into the calculated SAD [11]. Our 
results in Table 1 for the B c and B* meson masses are in a pretty good agreement with the 
other authers [1,4,7,11]. Thus, it is clear that the Song-Lin and Turin potentials are the 
most preferred interactions. Here, we report the range of the strong coupling constant at 
the m c scale we take in our analysis 0.1985 < a s (m 2 c ) < 0.320 for all types of potentials and 
0.220 < a s (ml) < 0.320 for the class of static potentials [11]. We point out a different choice 
of the potential would in general lead to a different value of the wave function at the origin 
and to a different determination of a s {m 2 c ) from the same hyperfine splitting. Further, our 
predictions to the be masses of the lowest S-wave (singlet and triplet) together with the other 
estimations by many authors are given in Table 2. Moreover, in Table 3, we also estimate 
the radial wave function of the low-lying state of the be system, so that 

|#is(0)| = 1.280 - 1.540 GeV 3/2 , (54) 

3 Kiselev et al. [34] have taken into account that AM T (1S) = ^4||AM V ,(1S) with a s (T) / a s (ip) ~ 
3/4. Further, Motyka and Zalewiski [21] also found ^ H/18. 



15 



for the group of static potentials. Further, we present our calculations for the NR leptonic 
constant f^f = 466^ MeV and /£, R = 463±24 MeV as an estimation of the potential 
models without the matching [4,19]. Our results are compared with those of other works 
[1,35,36]. Our calculation for the one-loop correction: 

^l-loop) = 40g +i6 ^l-loop) = 405 +l7 MeV; (55) 

and for two-loop correction: 

^2-loop) = 315 +16 ^2-loop) = 313 +26 Me y_ (5g) 

So, our numerical value for /]j R agrees with the estimates obtained in the framework of the 
lattice QCD result [5], /£ R = 440 ±20 MeV, QCD sum rules [37], potential models [1,4,19], 
and from the scaling relation [34]. It indicates that the one-loop matching [33] provides the 
magnitude of correction of nearly 12%. Further, the most recent calculation [29] in the heavy 
quark potential in the static limit of QCD with the one-loop matching is 

yO-ioop) = 40Q ± 15 MeV _ ( 57 ) 

Therefore, in contrast to the discussion given in [29], we see that the difference is not crucially 
large in our estimation to one-loop value in the 5 c -meson. On the other hand, our final result 
of the two-loop calculations is 

^2-loop) = 315 +26 MeV; (5g) 

the larger error value in (58) is due to the strongest running coupling constant in Cornell 
potential. Moreover, Motyka and Zalewiski [21] also found /s c ~ loop ' ) = 435 MeV for the 
ground state of be quarkonium. 

In the potential model, we note that slightly different additive constants is permitted 
to bring up data to its center-of-gravity value. However, with no additive constant to the 
Cornell potential [38], we notice that the smaller mass values for the composing quarks of 
the meson leads to a rise in the values of the potential parameters which in turn produces 
a notable lower value for the leptonic constant. 
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Our predictions for the be mass spectrum for the Igi-Ono potential (type I and II) are 
given in Table 4. Moreover, the singlet and triplet masses together with the hyperfine split- 
tings predicted for the two types of this potential are also reported in Table 5. We, hereby, 
tested acceptable parameters for from 100 to 500 MeV for the type I and II potentials 4 . 
Small discrepancies between our prediction and SAD experiment [11,16,19,23] can be seen 
for higher states and such discrepanciesare probably seen for any potential model and it 
might be related to the threshold effects or quark-gloun mixings. The fitted set of parame- 
ters for the Igi-Ono potential (type III) [11] are also tested in our method with b = 19 and 
(Ajfs = 300 MeV and also 390 MeV) and then b = 16.3 and = 300 MeV which seems 
to be more convenient than A^r = 500 MeV used by other authors [8]. Results of this study 
are also presented in Table 6. We see that the quark masses m c and m b are sensitive to the 
variation of A-^g. Therefore, as A^- increases the contribution of the potential (cf. e.g., 
Eqs. (31) and (32)) and consequently the binding energy E n j term decreases which leads to 
an increase in the constituent quark masses of the convenient meson, cf. Eq. (18). 

In this model, we see that the experimental be splittings can be reproduced for Ajjg ~ 
300 MeV in type I, A^ ~ 400 MeV in type II (cf. Table 5) and ~ 300 MeV in type 
III (cf. Table 6). We also predicted the splittings to several MeV with the other formalisms 
(cf., Table 1 of [6]). 

In Table 6, we also find that m c and nib are insensitive to the variation of Aj^ for this 
Chen-Kuang potential. This is consistent with the conventional idea that, for heavy quarks, 
the constituent quark mass is close to the current quark mass which is Ajjg independent. 
Numerical calculations show that this potential is insensitive to A^- in the range from 100 
to 300 MeV and as A.jjg increases, the potential becomes more sensetive for the lS-state 
only. The theoretically calculated n 1 So and n 3 Si hyperfine splittings for the B c meson in 
the Chen-Kuang potential are also listed in Table 6. They are considerably smaller than the 



4 The parameters of this potential are given in Table 3 of Ref. [11]. 
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corresponding calculated values A is (6c) = 76 MeV, and A 2S (5c) = 42 MeV predicted by 
the quadratic formalism of Ref. [6]. Moreover, Chen-Kuang [25] calculated A is (6c) = 49.9 
MeV, and A 2S (6c) = 29.4 MeV for their potential with = 200 MeV in which the last 
splitting is almost constant as increases. Our theoretical calculation for Ais(5c) = 68 
MeV, and A 2S (5c) = 35 MeV for the Chen-Kuang potential with A^ runs from 100 into 
375 MeV. We also find A is (6c) = 67 MeV, and A 2S (bc) = 33 MeV for the Igi-Ono potential 
with A^r = 300 MeV and b = 16.3. This model has the following features: (1) The 
present potential predicts smaller calculated A 1S and A 2 s than the other potentials do for 
be system and our calculated Ais and A 2 s do not depend on Ajjg more sensitively (2) The 
theoretical be splitting can be repoduced for the preferred A^ in the range 100 — 300 MeV. 
Furthermore, in Table 5, for instance, we may choose /|J R = 420 MeV with Aj^g = 300 MeV, 
for the type I, and /£ R = 396 MeV with A^ = 300 MeV, for the type II. The result on 
f Bc is within the errors given by the other authors [5,19,29]. Further, for the CK potential, 
we present the decay constant in Table 6. 

The scaling relation (SR) for the S-wave heavy quarkonia has the form [34] 

fl_ ( M n {bc) \ 2 ( m c + m b \ _ d 
M n {bc) \M 1 (bc)J { 4/x J n' 1 ' 

where m c and m& are the masses of heavy quarks composing the £? c -meson, fj, is the reduced 
mass of quarks, and d is a constant independent of both the quark flavors and the level 
number n. The value of d is determined by the splitting between the 2S and IS levels or 
the average kinetic energy of heavy quarks, which is independent of the quark flavors and 
n with the accuracy accepted. The accuracy depends on the heavy quark masses and it is 
discussed in detail [34]. The parameter value in Eq. (59), d ~ 55 MeV, can be extracted 
from the experimentally known leptonic constants of ip and T. So, from Table 1, the SR 
gives for the IS- level 

/£ R) ^ 444« MeV, (60) 

for all static potentials used. Kiselev [29,34] estimated f B = 400 ± 45 MeV and — 
385 ± 25 MeV, Narison [39] found /f c R) = 400 ± 25 MeV. 
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Overmore, we present the leptonic constants for the excited nS-levels of the be in Table 
7. We see that our prediction fg^s) = ± 15 MeV is in good agreement with the ones 
predicted by Kiselev et al. [19], /b c ( 2 s) = 280 ±50 MeV for the 2S-level in the be system. In 
Figure 1, we plot the calculated values of B c leptonic constants using Eq. (59) for different 
potential models together with the calculated values of the excited nS-states using [34] 

fn 2 = Jff ni . (61) 

The conclusion can be drawn from the Figure 1, that the approximated values of the excited 
nS-states agree well with the simple scaling relation (SR) derived from QCD sum rules for 
the state density. It is clear that the estimates obtained from the potential model and SR 
are in good agreement to several MeV as in Figure 1. Therefore, the difference between the 
leptonic constants for the pseudoscalar and vector lS-states is caused by the spin-dependent 
corrections, which are small. Numerically, we get f B » — f Bc / f B * < 1%. For the heavy 
quarkonia, the QCD sum rule approximation, provides that the fp and fy values for the 
pseudoscalar and vector states. leptonic constant is practically independent of the total spin 
of quarks, so that 

fv,n ^ fp,n = fn- (62) 

Our numerical approximation for the decay constants of the pseudoscalar and vector states 
in Tables 3, 5, and 6 is a confirmation to the last formula (62). 

In this paper, we have developed the SLET in the treatment of the be system using a 
wide class of static and QCD-motivated potentials. For such quarkonium potentials the 
method simply predicts the results of [11]. In this context, in reproducing the SAD, we used 
the same fitted parameters of [11] for the sake of comparison. Further, we demonstrate to 
the reader that the SLET method generates exactly the same numerical energy spectrum 
for cc, 66,and cb sresults as in the SLNET. This refutes the claims of [13] that this method 
is a reformation to SLNET and has a wider domain of applicability. Clearly, the method is 
simply an alternative parallel mathematical pseudoperturbative expansion technique having 
the same accuracy of [11,12,14]. 
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Appendix A: SLET Parameters for the Schrodinger Equation: 



Here, we list the analytic expressions of 7^ , 7^ , £j and 5j for the Schrodinger equation: 



7 M = 



(1 + 2n r )e 2 + 3(1 + 2n r + 2n 2 r )e 4 



LU 



ef + 6(1 + 2n r ) + (11 + 30n r + 30n;)e 2 3 



(Al) 



7 (2) = [(1 + 2n r )5 2 + 3(1 + 2n r + 2n 2 r )5 4 + 5(3 + 8n r + 6n 2 r + An 3 )5 & 

- LU' 1 (1 + 2n r ) e\ + 12(1 + 2n r + 2n 2 r )e 2 e 4 + 2e 1 5 1 
+ 2(21 + 59n r + 51n 2 r + 34n 3 r )e 2 4 + 6(1 + 2n r )e 1 S 3 

+ 30 (1 + 2n r + 2n 2 r ) ^65 +2 (11 + 30n r + 30n 2 r ) e 3 5 3 

+ 10(13 + 40n r + A2n 2 r + 28n 3 r )e 3 5 b + 6(1 + 2n r )£ 3 6i" 
+ uj' 2 he 2 e 2 + 36(1 + 2n r )e 1 e 2 e 3 + 8(11 + 30n r + 30n 2 ) e 2 e\ 
+ 24(1 + 2n r )e\e 4 + 8(31 + 78n r + 78n 2 )e 1 e 3 e 4 
+ 12 (57 + 189n r + 225n 2 r + 150^) e\ e 4 

- to' 3 \8ele 3 + 108(1 + 2n T )e\e\ + 48(11 + 30n r + 3^n 2 ^e x e\ 
+ 30(31 + 109n r + Uln 2 + Mn 3 r ) 4} , 



where 



(4/1^/2 ' 



% = 1,2,3,4. 



(A2) 



(A3) 



and 
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% = J^2> J = 1,2,3,4,5 ,6. (A4) 
_ ~(2a +1) _ 3(2a+l) 

£l " 2/x ' 52 _— i^ - ' (A5) 

63 = + "6^' £A = IH + 24Q (A6) 

ft = -^±±2; 5 2 = ^±1), (A7) 

* 3 = -^^; ^ 4 = ^^, (A8) 

3 7 r*V»»»(r ) 

65 = ~2jl + 120Q 5 5e= 4^ + 720Q " (A9) 
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FIGURES 

Fig. 1. The nS-levels leptonic constant of the B c system calculated in different static potential 
models using SR. 
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TABLES 

Table 1. The be masses and hyperfine splittings (A n s) calculated in different static potentials (in 
MeV). 



States 


Refs. 1,6 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


a s (m 2 c ) = 




0.320 


0.263 


0.286 


0.251 


0.220 


m c (GeV) = 




1.840 


1.820 


1.790 


1.800 


1.500 


m b {GeV) = 




5.232 


5.199 


5.171 


5.174 


4.905 


M(bc) 














IS 


6315 


6315 


6306 


6307 


6301 


6317 


1 3 S 1 


6334 


6335 


6325 


6326 


6319 


6334 


I'So 


6258 


6252 


6249 


6249 


6247 


6266 


A ls a 


77 


83.5 


76.1 


76.7 


71.6 


68.0 


2S 


6873 


6888 


6875 


6880 


6892 


6903 


2 3 S 1 


6883 


6897 


6884 


6889 


6902 


6911 


2 1 S 


6841 


6860 


6850 


6852 


6865 


6879 


A 25 


42 


37.9 


34.0 


36.5 


36.7 


31.3 


3S 


7246 


7271 


7209 


7246 


7236 


7225 


4S 




7587 


7455 


7535 


7483 


7448 


IP 


6772 


6743 


6733 


6731 


6730 


6754 


2P 


7154 


7138 


7104 


7123 


7125 


7127 


3P 




7464 


7371 


7428 


7398 


7375 


ID 


7043 


7003 


6998 


6998 


7011 


7027 


2D 


7367 


7340 


7284 


7320 


7311 


7301 


3D 




7636 


7510 


7588 


7536 


7502 



A nS = M(n 3 5i) - M^Sq). 
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Table 2. The predicted be masses of the lowest S-wave and its theoretically calculated splittings 
compared with the other authors (in MeV). 



Work 


M Bc (l l S ) a 


M B ,(l 3 5i) 




Eichten et al. [1] 


6258 ± 20 






Colangelo and Fazio [3] 


6280 


6350 




Baker et al. [40] 


6287 


6372 




Roncaglia et al. [40] 




6320 ± 10 




Godfrey et al. [1] 


6270 


6340 




Bagan et al. [1,40] 


6255±20 


6330 ± 20 




Brambilla et al. [3] 




6326tg 9 


60 6 


Baldicchi et al. [6] 


6194 ~ 6292 


6284 ~ 6357 


65 < Ais < 90 


SLET C 


6253i<t 3 


6328±^ 


68 < Ais < 83 


SLET d 


62581^ 


6333±? 4 





a The experimental mass of the singlet state is given in [2]. 

6 We remark that we have calculated splittings from the singlet and triplet states. 
c Averaging over the five values in Table 1. 

d We treat Eichten and Quigg's results in the same manner (e.g., [1]). 
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Table 3. The characteristics of the radial wave function at the origin |i?is(0)| 2 (in GeV 3 ), NR, 
one-loop and two-loop corrections to pseudoscalar and vector decay constants of the low-lying B c 
meson (the accuracy is 5%) alculated in different static potential models (in MeV). 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


GKLT[35] 


EFG[1] 


JW[36] 


l<M0)| 2 

\Rism 2 


0.112 
1.413 


0.123 
1.54 


0.111 
1.397 


0.119 
1.495 


0.102 
1.28 








ANR) 
J B c 

ANR) 

Jb* 


464.5 
461.5 


485.1 
482.2 


462.0 
459.2 


478.0 
475.3 


441.7 
439.3 


460±60 
460±60 


433 
503 


420±13 



loop) 

J B c 


393.6 a 


424.4 


399.6 


421.2 


399.3 


„(2~loop) 
J B c 


264.1 6 


333.0 


296.6 


339.1 


340.9 


j,(l—loop) 

Jb* 


391.0 


421.9 


397.1 


418.8 


397.2 


r-(2 — loop) 


262.3 


331.0 


294.8 


337.2 


339.0 



"First loop SD Wilson coefficient for all potentials, Kq = 0.85 — 0.90. 
^Second loop SD Wilson coefficient for all potentials, Kq = 0.57 — 0.77. 
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Table 4. The be mass spectra predicted for various Aj^g using Igi-Ono (type I and II) potential (in 
MeV). 

















States 


[6,24] 


100 


200 


300 


400 


500 


b = 20 a 


a s = 


0.1985 


0.217 


0.238 


0.250 


0.262 


15 


6327 


6329 


6318 


6310 


6316 


6327 


25 


6906 


6915 


6904 


6881 


6880 


6901 


35 


7246 


7264 


7242 


7244 


7241 


7252 


45 




7508 


7522 


7545 


7542 


7552 


IP 


6754 


6755 


6744 


6733 


6732 


6742 


2P 


7154 


7144 


7131 


7125 


7122 


7134 


ID 


7028 


7029 


7017 


7004 


7000 


7010 


2D 


7367 


7334 


7327 


7327 


7323 


7333 


b = 5 b 


a s = 


0.1985 


0.227 


0.230 


0.2405 




15 


6327 


6331 


6324 


6316 


6307 




25 


6906 


6914 


6898 


6910 c 


6918 




35 


7246 


7258 


7277 


7236 


7201 c 




45 




7521 


7517 


7478 


7500 




IP 


6754 


6756 


6743 


6737 


6730 




2P 


7154 


7142 


7138 


7134 


7120 




ID 


7028 


7029 


7015 


7012 c 


7007 




2D 


7367 


7335 


7323 c 


7314 


7316 





c = -0.022 to -0.031 MeV. 



b CQ = -0.019 to -0.026 MeV. 

c Carried out to the second correction order. 
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Table 5. The be mass spectrum, splittings and leptonic constant predicted for various Aj^g using 
Igi-Ono (type I and II) potential (in MeV). 



MS 



otates 


100 


200 


300 


A nn 
400 


500 


Type I 












l 3 5i 


6343 


6334 


6327 


6334 


6344 




6287 


6272 


6259 


6263 


6274 


Ais 


56.3 


62.0 


68.3 


71.1 


69.8 


\Ris(o)\ 2 


0.826 


1.005 


1.156 


1.19 


1.114 


fNR 
JB C 


9.^4 1 
oo^±. 1 


3Q1 1 


49D n 


49fi n 


41 1 7 


A\—loop) 
J B c 


one 1 


OOD.O 


oio.o 


o i y.z 


OD4.4 


r(2 — loop) 

JB C 


zyu.u 


one -i 

oUb.l 


oil 7 

oil. / 


onf! K 

oUb.o 


OC7 1 
ZO I .1 


fNR 

Jb* 


OOZ.O 


ooy.z 


A 1 7 7 
41 / . / 


4Zo.O 


/ino /i 


,.(1— loop) 

Jb* 


39fi 7 
OZO. / 


O04.O 


37/1 7 
O 1 4. / 


377 1 


3fi9 /I 
ODZ.4 


f(2— loop) 

Jb* 


288.7 


304.6 


o 1 n n 

310.0 


on yl T 

304.7 


285.0 


Type II 












l 3 Si 


6345 


6340 


6331 


6323 




1% 


6288 


6279 


6269 


6259 




Ais 


56.7 


60.6 


61.8 


64.4 




l^is(0)| 


0.819 


0.891 


1.03 


1.204 




fNR 

JB C 


352.7 


368.2 


396.0 


428.6 




Jl—loop) 

Jb c 


327.1 


334.9 


357.4 


382.1 




r(2 — loop) 


289.1 


283.0 


300.1 


314.4 




fNR 

Jb* 


351.2 


366.4 


394.1 


426.4 




r(l — loOp) 

Jb* 


325.6 


333.3 


355.7 


380.1 




j,(2—loop) 

Jb* 


287.8 


281.7 


298.6 


312.8 
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Table 6. The be mass spectrum, splittings and leptonic constant predicted for various Aj^g using 
Igi-Ono (type III) and Chen-Kuang potentials (in MeV). 



State 




IO (III) 






CK 




b = 


16.3 


19 


19 


5.1 


5.1 


5.1 




300 


300 


390 


100 - 300 


350 


375 


a s = 


0.250 


0.2505 


0.2205 


0.270 


0.270 


0.270 


IS 


6309 


6309 


6297 


6324 


6372 


6354 


2S 


6880 


6870 


6877 


6880 


6880 


6880 


3S 


7247 


7236 


7254 


7258 


6258 


6258 


4S 


7553 


7541 


7563 


7570 


7570 


7570 


IP 


6725 


6721 


6737 


6723 


6723 


6723 


2P 


7124 


7114 


7135 


7127 


7127 


7127 


3P 


7441 


7429 


7452 


7452 


7452 


7452 


ID 


6997 


6990 


7013 


6993 


6993 


6993 


2D 


7328 


7317 


7341 


7332 


7332 


7332 


3D 


7613 


7599 


7624 


7625 


7625 




1 3 5! 


6326 


6327 


6315 


6341 


6389 


6371 


I 1 S 


6259 


6258 


6243 


6273 


6321 


6304 


Ais 


67.3 


68.6 


72.6 


67.8 


67.8 


67.7 


I^15(0)| 2 


1.115 


1.119 


1.339 


1.017 


1.017 


1.017 


fNR 
JB C 


412.4 


413.2 


452.6 


393.5 


392.0 


392.5 


A\—loop) 
J B c 


367.1 


367.2 


408.0 


347.5 


346.2 


346.6 


r(2 — loop) 

JB C 


296.4 


294.0 


345.4 


269.0 


268.0 


268.4 


fNR 

Jb* 


410.2 


411.0 


450.0 


391.4 


390.0 


390.5 


r(l — loOp) 

Jb* 


365.2 


365.2 


405.7 


345.6 


344.3 


344.8 


j,(2—loop) 


294.8 


292.4 


343.4 


267.6 


266.6 


266.9 
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Table 7. The nS-levels leptonic constant of the be system, calculated in different static potential 
models (the accuracy is 3 — 7%), in MeV, using the SR. 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


/is 


449.6 


450.4 


448.0 


448.8 


420.9 


/2s 


305.8 


305.0 


303.3 


303.5 


284.7 


/3s 


243.0 


243.2 


241.3 


241.8 


227.2 


/4s 


206.0 


207.1 


204.9 


205.9 


193.8 
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